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Testing dark energy using pairs of galaxies in redshift space 
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ABSTRACT 

The distribution of angles subtended between pairs of galaxies and the line of sight, 
which is uniform in real space, is distorted by their peculiar motions, and has been 
proposed as a probe of cosmic expansion. We test this idea using N-body simulations of 
structure formation in a cold dark matter universe with a cosmological constant and in 
two variant cosmologies with different dark energy models. We find that the distortion 
of the distribution of angles is sensitive to the nature of dark energy. However, for 
the first time, our simulations also reveal dependences of the normalization of the 
distribution on both redshift and cosmology that have been neglected in previous 
work. This introduces systematics that severely limit the usefulness of the original 
method. Guided by our simulations, we devise a new, improved test of the nature of 
dark energy. We demonstrate that this test does not require prior knowledge of the 
background cosmology and that it can even distinguish between models that have 
the same baryonic acoustic oscillations and dark matter halo mass functions. Our 
technique could be applied to the completed BOSS galaxy redshift survey to constrain 
the expansion history of the Universe to better than 2%. The method will also produce 
different signals for dark energy and modified gravity cosmologies even when they have 
identical expansion histories, through the different peculiar velocities induced in these 
cases. 

Key words: Methods: N-body simulations - Cosmology: theory - large-scale structure 
of the Universe - dark energy 



1 INTRODUCTION 



One of the primary scientific goals of ongoing and future 
galaxy surveys is to determine what i s responsible for the 
accelerating expans i on of t he Universe jLaureiis et al.l 2010l: 
Schlegel et all 120071. 120091 : iBlake et al.l l2010l : ISchlegel et al l 



201lft . There are two main considerations which affect cur- 



rent tests of dark energy. Firstly, the differences between the 
observables expected from competitive cosmological models 
are small. Secondly, given the huge volumes that will be 
covered by future surveys it is likely that systematic errors 
will dominate the interpretation of the measurements. For 
both these reasons it is generally accepted that the dark en- 
ergy ch allenge should be tackled using multiple cos mological 
probes l|Peacock et al.ll2006l : lAlbrecht et all 12006ft . Guided 
by numerical simulations, we assess a re cently proposed test, 
a ver sion of the Alcock-Paczynski test (|Alcock fc Paczvnskil 
which uses measurements of galaxy pairs to constrain 
the cosmological model. We expand the available probes of 
dark energy by setting out an improved version of the test 
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which we show can distinguish models that otherwise cannot 
be separated by existing methods. 

The Alcock-Paczynski test measures the distortion of a 
spherical object assuming an incorrect cosmological model 
is used to compute distances. The version of this test con- 
sidered here models the distortion in a spherically symmet- 
ric distribution of galaxy pair angles in redshift space and 
was first proposed in a form similar to that used in this 
paper by iPhillippsj i| 1994ft . who considered the distribution 
of ang les between quasar pairs. Recently iMarinoni fc Buzzil 
l|2010ft introduced an important revision to this test by con- 
sidering the angle between pairs of galaxies viewed in red- 
shift space. Building on the work of Marinoni & Buzzi, the 
method outlined in this paper is a refined geometrical test 
of dark energy. The critical feature of our extension is the 
use of N-body simulations of different dark energy models to 
test the idea that measuring the anisotropic distribution of 
galaxy pairs in redshift space is a useful probe of cosmology. 
This geometrical test of dark energy will complement and 
extend currently used geometrical probes such as measur- 
ing the light curves of Type l a supernovae (SN) (|Riess et al.l 
1 19981 : iPerlmutter et alj ll999l ; iRiess et all l201lft and appli- 
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cations of the Alcock-Paczynski test to baryonic acoustic 
oscillations (e.g. Blake et a l. 201 it ) ■ 

In a Friedmann-Robertson- Walker universe, pairs of 
galaxies should be distributed with random orientations if 
the fundamental assumptions of homogeneity and isotropy 
are correct. This simple test of cosmology is complicated by 
two effects: firstly, we do not observe galaxies directly in real 
space but in redshift space, where peculiar velocities, dis- 
tinct from the Hubble flow, displace the position of a galaxy 
along the line of sight from its true position. This introduces 
a preferred direction, with the result that galaxy pairs are 
no longer randomly distributed. Secondly, in order to con- 
vert observed angles and redshifts into comoving distances, 
an observer needs to assume a cosmological model. An in- 
trinsically spherical object, such as a cluster of galaxies, or 
a spherically symmetrical distribution, such as the distri- 
bution of galaxy pairs we consider in this paper, will ap- 
pear distorted if measured assuming a cosmology that does 
not match the true under lying cosmology of the Universe 
l|Alcock fc Paczvnskj|l979h . 

Based on this idea, IPhillippsI (|l994T ) proposed 
a test where the hyp othetical sphere proposed by 
lAlcock fc Paczvnskil (|l979l ) is replaced by randomly orien- 
tated quasar pairs. In the absence of peculiar motions, a 
large sample of quasar pairs should have a uniform distri- 
bution in the cosine of the a ngle between each pair, if the 
correct cosmology is adopted. iMarinonf fc Buzzil i|2010l ) de- 
veloped this test by modelling the effect of the redshift space 
distortions as a Doppler shift in the positions of the galaxies. 
They applied this mo del to galaxy pairs in t he Sloan Digital 
Sky Survey (SDSS) (|Abazaiian et al.ll2009l) at low reds hift 
and the DEEP2 galaxy redshift survey (|Coil et al.ll2"004h at 
z ~ 1.3. After selecting galaxy pairs according to a set of 
constraints discussed in Section 4, Marinoni & Buzzi were 
left with a sample of 721 pairs in the SDSS DR7 at z ~ 0. 
Mainly due to the small sample size, Marinoni & Buzzi were 
only able to distinguish a ACDM cosmology from somewhat 
extreme alternatives, namely an Einstein-de Sitter Universe 
and an open universe with no dark energy. We note that 
these models have already been ruled out by other tests. 
This does not, however, imply that the test cannot be used 
to yield competitive constraints on dark energy with a larger 
sample of pairs. 

Note that the Alcock-Paczynski test measures a distor- 
tion parameter which is proportional to the angular diame- 
ter distance at the redshift of the object, Da(z) multiplied 
by the Hubble rate, H(z). The test we propose in this paper 
models the distribution of galaxy pairs in real and redshift 
space assuming a distant observer approximation. This as- 
sumption removes any dependence on the angular diameter 
distance in the distortion parameter. As a result this tech- 
nique allows us to measure the Hubble rate directly. 

In this paper, using subhalo pairs in large volume N- 
body simulations, we test the method of Marinoni & Buzzi 
and its potential to distinguish between cosmologies. First 
we focus on the selection criteria necessary to provide a ho- 
mogeneous sample of pairs whose distribution in redshift 
space agrees with the theoretical model of Marinoni & Buzzi. 
With robust selection criteria, we then apply this test to dif- 
ferent simulations to see if these dark energy models can be 
distinguished from ACDM. A critical assumption made in 
the analysis by Marinoni & Buzzi is that the normalization 



of a theoretical model of the pair distribution does not evolve 
with redshift. We show, using numerical simulations, that 
this assumption is incorrect. We also consider the practical 
difficulties associated with obtaining an accurate measure- 
ment of this normalization parameter observationally. We 
demonstrate that the test, as originally proposed, suffers 
from large systematics which limit its utility. We present an 
improved methodology which uses N-body simulations and 
does not require prior knowledge of the true cosmological 
model. 

Most attention to date has focused on cosmological tests 
which require measurements on large scales, such as the 
rate at wh ich cosmic structures grow (|Guzzo et alj 120081 ; 
Wang 200 8]), the apparent loca t ion of baryonic ac oustic os- 
cillations (|Sanchez et all 120091 ; iBlake et all |2010| ) and the 
projected matter de nsity as measured through weak lensing 
jMassev et al.l l2007). It is important to expand this arsenal 
of tests. This introduces sensitivity to different systematics, 
which, alongside results from other probes, will lead ulti- 
mately to a more convincing measurement of the properties 
of dark energy. Also, it is useful to devise new tests which 
are not reliant on measuring the galaxy distribution on the 
very largest scales, thereby avoiding the need for an accu- 
rate determi nation of the mean galaxy density (for another 
example see iNusser et al.l l201lh . The test proposed in this 
paper requires a large volume simply to accumulate a large 
sample of galaxy pairs; there is no requirement implied on 
the accuracy of the photometry across a survey used for this 
purpose. 

A novel feature of our analysis is the use of N-body 
simulations to validate and improve upon the method- 
ology proposed by Marinoni & Buzzi. Recent work has 
shown that numerical simulations of structure formation 
have an important role to play in modelling cosmologi- 
cal probes and assessing potential systematic errors. An- 
gulo et al. (2008) demonstrated that the shape of the 
power spectrum of galaxy clustering is substantially dif- 
ferent from the predictions of lin ear perturbation t heory 
even on very la r ge scales (see also ISmith et al.l 120071 . 120081 ; 
ISeo et aTl 120081 ; I Jennings et al.l l2010t ). The simulation re- 
sults led to revised analyses of baryonic acoustic oscillations, 
which either attempt to model the distortions introduced 
by nonlinearities an d redshift space, or to reconstruct the 
linea r theory signal dSeo Sz Eisens tcin 2007; I Sanchez et al.l 
|200S| ; iMontesano et al.ll2010l . l201ll ). Similarly. N-body sim- 
ulations have demonstrated that the measurement of the 
growth fact or from redshift space distortions requires carefu l 
modelling djennings et all l2011al ; lok umura fc Jing|[20TTI l. 
I Jennings et all (|2011bl f showed that a naive application of a 
linear theory model for the distortion of clustering in red- 
shift space can lead to a catastrophic misinterpretation of 
the measured growth factor. The study in this paper is in 
a similar spirit; the availability of N-body simulations to 
model the pair distribution allows us to devise an improved 
cosmological probe. 

The outline of this paper i s as follows: In Section 2 , we 
review the theoretical model o f M armom fc Buzzil (|2010h for 
the anisotropic distribution of pairs and its dependence on 
cosmology. In Section 3 we discuss the quintessence dark en- 
ergy models and the N-body simulations used in this paper. 
In Section 4 we list and test the selection criteria used to se- 
lect a homogeneous sample of galaxy pairs from the N-body 
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simulations whose distribution agrees with the theoretical 
predictions. In Section 5 we present our results, comparing 
the theory with measurements from simulations, for the two 
dark energy and the ACDM cosmologies and demonstrate 
that a robust test of cosmology can only be achieved by 
combining observations with numerical simulations. In Sec- 
tion 6 and 7 we present our summary and conclusions. 



THEORETICAL BACKGROUND: THE 
DISTRIBUTION OF GALAXY PAIRS 



Follow ing the derivation and discussion in lMa rinoni & Buzzil 
(2010) (see also the alternative derivation in |Phillippd[l994h . 
let us consider a pair of gravitationally bound galaxies, A 
and B, at an observed angular separation, 8, as shown in 
Fig. Q] In a flat universe, the angle that galaxy B subtends 
at galaxy A, as measured from the observer's line of sight 
through A, which we refer to as the tilt angle, t, can be 
written as 



= [1 + cot 9 



XA 



Xb sin J 

(see lMarinoni fc Buzzil l|2010h for the general expression in a 
curved universe). Here Xa(b) is the radial comoving distance 
to galaxy A(B) which is given by 



dz' 



aoJo H(z') 



(2) 



where a = 1/(1 + z) is the scale factor of the universe with 
current value do = 1, c is the speed of light and H — a/a 
is the Hub ble parameter with current value Hq = 71.5 
km/s/Mpc (Sa nchez et alJl2009l ). In an isotropic and homo- 
geneous universe, the orientation of pairs of bound galax- 
ies will be randomly distributed, allowing us to predict the 
probability distribution for t, F(t), and a measure of the dis- 
tribution /i = (s'm 2 t). In the absence of peculiar velocities, 
an observer calculating the ensemble average of Eq.[T]should 
find a value of fi = 2/3, independent of cosmology, as long 
as the correct cosmology is assumed when converting angles 
and redshifts to comoving distances for each member of the 
pair. 

Galaxies have velocities in addition to the Hubble flow 
which result in inferred positions which appear displaced 
along the line of sight. As a result the true angle t between 
a pair of gravitationally bound galaxies will appear as an 
angle, r (see Fig.[T]). The angle subtended between the pair 
of galaxies in redshift space calculated according to Eq. [T] 
corresponds to a measurement of the angle r and not t, and 
the average over all pairs will no longer be a random distri- 
bution. The result is a skewed distribution whose mean will 
differ from the expected value for a uniform distribution of 
H — 2/3. Adopting the correct cosmology to calculate the 
ensemble average in Eq. Q] will then provide a measure of 
the mean of the apparent distribution in redshift space af- 
ter intrinsic peculiar velocities distort the orientation of the 
pairs. Marinoni & Buzzi modelled this distortion as a simple 
Doppler shift where the observed line of sight separation is 
related to the actual separation, to first order in v/c, by 



dr D 



dr 



dv \\ 
H(z) 



(l + z) : 



(3) 




Figure 1. An illustration of the actual, t and observed, t, angles 
subtended between a pair of galaxies, A and B, and the line of 
sight through A in real space (black) and in redshift space (red) 
respectively. At the observer's position, the pair subtend an angle 
8. The comoving line of sight separation of the pair is shown as 
Arii and Arii obs in real and redshift space respectively. 



where dr oba and dr are the observed and actual line of 
sight separations of a pair of galaxies A and B and dv\\ = 
va ■ ta — vb ■ tb is the line of sight peculiar velocity dif- 
ference, where ta(b) represent unit vectors in the direction 
of each galaxy in the pair. The line of sight comoving sep- 
aration is dr = dx, as given in Eq. [2] Note this equation is 
a result of relating the position observed in redshift space 
to the actua l posit ion in real space as dz ohs = dz + dv\\ 
e.g lHamiltonl (|l998l ). In the distant observer approximation, 
the separation between galaxies is assumed to be small com- 
pared to the distance between them and the observer. Under 
this assumption the observed comoving transverse separa- 
tion is equal to the true transverse separation of the pair, 
drx i0bs ~ dr±. If the redshift difference of the pair Az is a 
lot less than unity such that Az ~ dz, then the observed, r, 
and actual tilt, t, of the pair can be simply related by the ob- 
served, Arii obs and the actual line of sight finite separation, 
Arn , according to 



tan t 
tanr 



Arii 



Am 



1 + 



(1 + z) Ai-n 
H(z) Ary 



(4) 



The relation given in Eq. [4] can then be used to transform 
the true distribution of galaxy pairs, F{i), into the apparent 
distribution, 9(t). Using conservation of probability, Mari- 
noni & Buzzi derived the probability distribution function 
of the apparent angle written in terms of the true angle as 



*(r)dr = F(t)dt. 
From this it follows that ^(r) is given by 
1 (l + cr 2 )(l + tan 2 r 



*(r)dr 



2 (1 + (1 + cr 2 ) tan 2 r 



.3/2 



tan r Idr . 



(•») 



(6) 



and the parameter a depends on the cosmological expansion 
history as 



.2... ^ o/ A ^ll\ l + « , ^Hi(l + z) 2 



a\z,n) = 2(-^ 



\ Ar I H(z) 



(7) 
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The normalization parameter a is given by 

2 \ \ V2 



= Hn 



Ar 2 



(8) 



The first moment of the distribution ^(t), referred to by 
Marinoni & Buzzi as the 'average anisotropy of pairs' ( the 
AAP function from now on), is given by 



Mobs 



(1 + a 2 ) arctan(cr) — a 



(9) 



In Eq. [9l the parameter a depends on the expansion history 
in a particular cosmological model, H(z), as given in Eq. 
\7\ Marinoni & Buzzi set the first term on the right hand 
side of Eq. [7] to zero on the assumption that the comoving 
separation of pairs and their radial peculiar velocities are un- 
correlated, (A«ii/Ar) = 0. We shall discuss this assumption 
further in Section 4. 

The original Alcock-Paczynski test, when applied to a 
spherical object, measures a distortion parameter, the ratio 
of the tangential and radial distances, which is proportional 
to Da(z)H(z) and is unity if the correct cosmological model 
is assumed and there are no redshift space effects i.e. there is 
no distortion of the spherical object. The Alcock-Paczynski 
test applied in this paper compares the distribution of pair 
angles in real and redshift space in the distant observer ap- 
proximation, drx, obs ~ drx, which gives rise to a distortion 
parameter which is independent of Da- The distortion is es- 
timated, after modelling redshift space effects, by comparing 
the distribution of the angles t and r, and only depends on 
H(z). 

Using pairs of galaxies in a survey, an observer can 
measure the average orientation using Eq. [1] which should 
be equal to the AAP function in Eq. [9] if the correct cos- 
mology is assumed and the observer is able to measure a 
precisely in order to fully specify ^(t). In this paper we 
perform this exact test using pairs of subhaloes in N-body 
simulations of different cosmologies. In practise in a galaxy 
survey the parameter a can be determined in two ways: 
firstly, at low redshifts, where the peculiar velocities of the 
pair can be measured by combining a redshift independent 
distance measurement, e.g. luminosity distances from Type 
la supernovae, the Tully-Fisher relat i on or the D n — cr re- 
lation (see e.g. ICourteau et al.l |2000|; Ida Costa et all l200d : 



iBorgani et aitl200d ; iBlakeslee et alj l200Ch . with the mea- 
sured redshift of the galaxy. The uncertainties associated 
with the redshift independent luminosity distance measure- 
ments are large, ~ 10 — 20% for the Tully-Fisher or D„ — a 
relations, and ~ 5—10% f or supernovae or the surface bright- 
ness fluctuation method l|Bernardi et al. I l200ll ; iTonrv" et al.l 
2001). These uncertainties on the redshift independent dis- 
tance measurements propagate into larger errors for the pe- 
culiar velocities, making it almost impossible to accurately 
measure the peculiar velocity of a single galaxy. The sec- 
ond method to determine a observationally, which we shall 
assess in this paper after considering the ideal case of mea- 
suring a from the simulations using Eq. 8, is to fit to the 
measured distribution of pairs at each redshift using Eq. [S] 
One of the key assumptions made by Marinoni & Buzzi 
is that the normalization factor a is constant for all red- 
shifts and for different galaxy selections. At z ~ 0, Mari- 
noni fe Buzzi obtained a = 5 .791"°^, using binaries in the 
SDSS (|Abazaiian et all 120091) . Marinoni & Buzzi obtained 



this value by fitting Eq. [6] to the observed distribution at 
2 ~ 0. We explicitly test this assumption in this paper where 
it is possible to measure a directly from the N-body simula- 
tions at each redshift. We can also compare the predictions 
of the AAP function using the best fit value for a obtained 
at z = 0, instead of normalizing the function at each red- 
shift. This will allow us to see if the value of a really is 
independent of redshift. 



3 TRIAL SAMPLES OF PAIRS FROM 
NUMERICAL SIMULATIONS 

As a test of the method proposed by Marinoni & Buzzi which 
was outlined in Section 2, we apply it to different cosmolo- 
gies, focussing on quintessence models. In Section 13.11 we 
discuss the two quintessence dark energy models we take as 
examples and highlight the main differences between these 
and the concordance cosmological model. In Section 13.21 we 
describe the simulations carried out. 



3.1 Quintessence dark energy 

Numerous quintessence dark energy models have been 
considered a s an alternative to the concordance cosmol- 
ogy (see e.g. Ratra fc Peebleslll988l ; iFerreira fc Joyce! Fl998l ; 
ICopeland et al.ll200d ; lMartmll2008f ). We focus on two inter- 
esting examples which are representative of a wider class 
of quintessence models, scalar fields which evolve in time, 
which are viable alternative cosmologies. 

One of the models we consider has substantial differ- 
ences to ACDM and can be considered as an 'early' dark 
energy model which features non-negligible amounts of dark 
energy at high redshifts. This quintessence dark energy 
model features an exponential term in the scalar field po- 
tential which pushe s the dark energy eq uation of state to 
w = -0.82 today (|Brax fc Martini 1 19991 ) . We refer to this 
model as the SUGRA model. The second quintessence dark 
energy model, which we refer to as INV, has been shown to 
produce a similar expansion history and no n-linear growth 
of structure to those in a ACDM cosmology (|jennings et al.l 
I20ld . l2011ah and will provide a measure of the sensitivity of 
the test we perform in this paper. The INV model has an in- 
ve rse power law poten tial V{<j>) = A l3+i /(j> for the scalar field 
(f> (|Zlatev et al.l 1999). The values of the constants A and /3 
are fixed by the current value of t he dark energy density (see 
e.g. ICorasaniti fc Copelandll2003l ). 

The dark energy equation of state for these quintessence 
models can be accurately des cribed over a wide range of red - 
shifts using four parameters l|Corasaniti fc Copelandl feoOS'). 
The variables used are: wq, the current dark energy equation 



of state; w m , the value of w during the matter dominated 
the scale factor at which the dark energy equation 
of state changes from its value during the matter dominated 
era to its present value, and A m , the width of the transi- 
tion in the expansion factor. For the SUGRA model these 
parameters are Wo = — 0.82, w m = — 0.18, a m = 0.1 and 
A m = 0.7. For the INV model the values of the parame- 
ters are wp 
l|jennings et al.ll20ld l 

The dark energy models have different expansion his- 
tories to ACDM and so when compared to the currently 



0.79, w m = -0.67, a m = 0.29 and A m = 0.4 
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ng_ 

ues of the cosmological parameters (see lJennings et al.ll20ld . 
for a discussion). As our starting point, we consider a 
ACDM model with the following cosmological parameters: 
f2 m = 0.26, Ode = 0.74, fi b = 0.044, h = 0.715, where 
Hp = lOOfe km/s/M pc and a spectral tilt of ns = 0.96 
jSanchez et al l [2009 ) . The linear theory rms fluctuation in 
spheres of radius 8 Mpc is set to be as — 0.8. In the 
simulations discussed in this paper, the ACDM values for 
fl m and Ho were used for the INV dark energy model while 
for the SUGRA model the best fit parame ters used were 
Qm = 0.243 and Ho = 67.73km/s/Mpc (see I Jennings et all 
2010, for more details). Both of these models are con- 
sistent with current observation s of Supernovae Type la 
lightcur ves (iKowalski et al.ll2008|). barvonic acou stic oscil- 



lations |Percival et al. 20071 ; Sanchez et al.l 120091 ') and the 



WMAP 7 year measurements of the cosmic microwave back- 
ground |Komatsu et al.l f2010). A detailed study of both of 
th ese models compared to a ACDM cosmology can be found 



in I Jennings et al" (l2010h. 

Note I Jennings et all |2010l ) showed that the INV model 
was indistinguishable from ACDM for several cosmological 
probes such as measurements of the halo mass function, 
BAO peak positions and growth factor. This model provides 
us with a significant test of the discriminatory power of the 
technique proposed in this paper. 

The ratio of the Hubble parameter in each quintessence 
dark energy model to that in a ACDM cosmology is shown 
as a function of redshift in Fig. [21 



3.2 N-body simulations 

The simulations were carried out at the Institute of Compu- 
tational Cosmology using a memory efficien t version of th e 
TreePM code Gadget-2, called L-Gadget-2 (|Springelll2005r) . 
The simulations use TV = 646 3 ~ 269 x 10 6 particles to rep- 
resent the dark matter in a computational box of comoving 
length 1500/i _1 Mpc. We shall refer to these simulations as 
the low resolution runs in Section [57T] We chose a comoving 
softening length of e = 50/i _1 kpc. The particle mass in the 
simulation is 9.02 x 10 h Mq with a mean interparticle 
separation of r ~ 2.3 ft _1 Mpc. We also consider a higher res- 
olution simulation of the ACDM cosmology with the same 
boxsize as above but with 1024 3 dark matter particles, ap- 
proximately four times the number of particles used in the 
lower resolution simulation. 

The initial conditions of the particle l oad were set 
up with a glass configuration of particles (|Whitel 1 19941 ; 
iBaugh et al.l 19951 ). The particles are perturbed from the 



glass using the Zel'dovich approximation which can induce 
small scale transients in the measured power spectrum. 
These transients die away after ^3 — 10 expansion factors 
from the starting redshift (|Baugh et al. I ll995l ; ISmith et al.l 
2003). In order to limit the effects of the initial displace- 
ment scheme we chose a starting redshift of z — 200. The 
linear theory power spectrum used to generate the ini- 
tial conditions was c reated using the CAMB package of 
iLewis fc Bridle] (|2002h . The linear theory power spectrum 
for the two dark energy models was compute d using the 
Para meterized Post-Friedmann (PPF) module (Fang et al. 
2008) for CAMB which takes into account the effects of a 
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Figure 2. The ratio of the Hubble rate in the two quintessence 
dark energy models simulated to that in a ACDM cosmology 
plotted as a function of redshift. Note the SUGRA quintessence 
model has a current value for the Hubble parameter of Hq = 
67.73km/s/Mpc, consistent with observations (see Section 3.1 for 
details). 



dynamical dar k energy equation of state and dark energy 
perturbations II Jennings et al.ll2010t ). 

Dark matter haloes were identified in the simulation 
outputs using a friends-of-friends (FOF) percolation algo- 
rithm with a linking le ngth of b = 0.2 t imes the mean in- 
terparticle separation_m^^£aD [l985|) ■ The SUBFIND 
algorithm (Spr ingel et al.ll200ll ) was then run on these halo 
catalogues to identify self bound subhaloes at each redshift. 
Note that the subhaloes are not necessarily bound to the 
main FOF halo. In this paper, pairs of subhaloes within a 
common FOF halo are used as a proxy for pairs of galaxies. 

The minimum number of particles per halo and subhalo 
is 10 and we select only haloes that have at least two sub- 
haloes (i.e. a minimum of 20 particles in the FOF group). 
In Fig. 6 we show that our results are not affected by our 
choice of minimum FOF halo mass. 

The position of each subhalo in redshift space is com- 
puted by perturbing its comoving position in real space us- 
ing the line of sight centre of mass velocity of the subhalo 
relative to an observer placed at the origin of the box. At a 
given redshift z > the observer is still assumed to be at 
the origin of the box at z = which requires us to add the 
comoving distance from z — to z to the subhalo positions. 
In Section 4 we discuss the selection criterion for subhalo 
pairs in relation to the radius -R200 of the parent halo where 
the mean density is 200 times the critical value. 



4 CALIBRATING THE METHOD USING 
SIMULATIONS 

From the N-body simulations we know which pairs of sub- 
haloes are in the same FOF halo. However, this does not 
guarantee that these objects are gravitationally bound to 
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Figure 4. Left panel: The distribution of the line of sight peculiar velocity difference, Any , for all pairs of subhaloes in the lower 
resolution ACDM simulation at z = 0. Pairs to the left of the grey shaded region represent those subhaloes that have been selected (95% 
of the distribution). Lower right panel: The comoving transverse separation in redshift space, Ar±, of pairs of subhaloes in the lower 
resolution ACDM simulation at z = plotted as a fraction of ffaoo- The selected subhaloes are shown as a red hashed region while those 
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region. 
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Figure 3. The distribution of the cosine of the apparent incli- 
nation angle, cos(r), of pairs of subhaloes in a ACDM simulation 
at z = 0. The distribution of cos(r) for all the subhalo pairs in 
redshift (real) space is shown in red (black). Subhalo pairs which 
are selected according to the criteria discussed in Section 4 are 
shown as blue hashed boxes. Selecting pairs with only a cut in 
#max gives rise to the distribution shown in orange. Selecting pairs 
with only a cut in Ai> max gives rise to the distribution shown in 
green. 



the FOF halo. There are too many haloes in our simulations 
to check explicitly for binding, so we will use proxies instead. 
This will allow us to make contact with the observational se- 
lection applied by Marinoni & Buzzi and to see how their 
cuts translate into cuts in simulation quantities. We inves- 
tigate these selection criteria and provide robust selection 
cuts which are independent of cosmology and redshift. 

In particular we address the following question: how do 
we construct a sample of pairs that matches the theoretical 
expectation for the AAP function, in the most favourable 
case in which we know the correct cosmology? Using infor- 
mation output from the simulations about the subhaloes se- 
lected and the properties of the parent halo, (e.g. the FOF 
algorithm returns R200), we can quantify the definition of 
a close pair in a rigorous way. If we selected only bound 
pairs we would expect good agreement with the AAP func- 
tion, provided that we know the correct cosmological model. 
As a first approach to identify suitable pairs, we shall se- 
lect subhaloes which are within -R200 of the main halo, i.e. 
Arx.max = R200, with no other restrictions on velocity or 
distance to a nearest neighbour. As this information is not 
available to an observer, our second approach will be to 
translate these selection criteria into observable quantities 
such as the angle 9 between a pair of subhaloes. 

Marinoni & Buzzi used the following selection criteria 
to pick their sample of pairs: (1) a maximum line of sight 
velocity difference of the pair AV = 700 km/s to avoid pro- 
jection of neighbouring systems, (2) a maximum comoving 
transverse separation of Arx, max = 0.7 Mpc//i, (3) a min- 
imum comoving transverse separation Arj_ imln = 20 kpc/h 
and (4) a minimum comoving distance from the centre of 
the galaxy pair to another galaxy. The latter two conditions 
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avoid selecting pairs which may be in the process of merg- 
ing or which are interacting with another galaxy. The value 
for the maximum velocity difference was chosen such that 
the relative increase AN/N in the sample size was < 1% 
when the velocity cut was increased by 100 km/s, while the 
maximum comoving transverse separation was chosen to be 
equal to the distance from Andromeda to the Milky Way. 

Fig. [3] shows the measured distributions of the orienta- 
tion of subhalo pairs in real and redshift space in the low 
resolution ACDM simulation at z = 0. The starting point 
is the sample of subhalo pairs within a common FOF halo, 
without any further selection. This is shown in real space by 
the black histogram in Fig. [3] Note that for the lower res- 
olution ACDM simulation there are approximately 65,000 
subhalo pairs at z = 0. The real space distribution of the 
tilt follows the expected random distribution and is uniform 
in cos(r). The distribution of all subhalo pairs in redshift 
space is shown in red which is clearly skewed. The mean of 
this distribution differs from the prediction of the AAP func- 
tion by ~ 40%. Applying the final set of cuts to this overall 
sample as outlined below leaves approximately 19,000 pairs, 
and produces the blue hashed region which is skewed to- 
wards smaller angles and agrees with the predictions of the 
AAP function given in Eq. [9] to within 0.5%. We discuss 
the selection cuts that give rise to this blue hashed region 
below. A comparison of the red and blue histograms in Fig. 
[3] shows that in redshift space if no selection cuts are made 
to isolate bound pairs the distribution is clearly randomized 
by outliers. 

In an attempt to isolate subhaloes that are gravitation- 
ally bound to their parent FOF halo and hence to test if 
their orientations in redshift space are distributed accord- 
ing to the predictions of the AAP function, we first select 
pairs of subhaloes within R200 and exclude all other pairs. 
We find that this sample of pairs has a non-negligible cor- 
relation between Awy and Ar such that (Av^/Ar) 7^ 0. As 
a result, we use the full expression in Eq. [7] for the param- 
eter a. This gives an AAP function which is in remarkably 
good agreement with the measured mean of the distribu- 
tion, the ensemble average of Eq. 1, at z = in a ACDM 
simulation, to better than a percent. This agreement dimin- 
ishes at higher redshifts, with the AAP function and the 
measured mean differing by 10-20% over the redshift range 
z = 0.25 - 1. 

It is possible to remove subhalo pairs which have 
(A«||/Ar) 7^ by selecting pairs according to an upper limit 
in the line of sight peculiar velocity difference, Ai> max . The 
velocity difference of pairs of galaxies is related to the com- 
mon gravitational potential of the pair which, in most cases, 
is weakly correlated with their separation. However, we find 
that pairs with large velocity differences have non-zero cor- 
relations, e.g. in the ACDM simulation at z = 0, using all 
subhalo pairs with Av > 950km/s we find (Av^/Ar) = 8.5h 
km/s/Mpc. Observationally these subhaloes would not be 
detected as the apparent tilt between the pair is approxi- 
mately zero, due to their large peculiar velocity difference, 
and as a result the pair would lie along the same line of sight. 
In Fig. [4] we plot the distribution of the line of sight pecu- 
liar velocity difference Ad in the left panel, for all subhaloes 
in the lower resolution ACDM simulation. The grey shaded 
region corresponds to the selection cut in Av. Once we re- 
move any correlated pairs from the sample, and impose the 



restriction that Arx jmax = -R200, we find that the measured 
mean and the predicted AAP function agree extremely well 
in the redshift range 2 = — 2. We present these results 
in more detail in the following section. Note the first term 
in the expression for a, Eq. [7] is now negligible as we have 
removed any correlated pairs. 

As -R200 is not an observable quantity, the next step is 
to see if this cut can be translated into a cut in 9, the ob- 
served angular separation of the pair. In Fig. U we plot the 
distribution of the comoving transverse separation Ar± as a 
fraction of -R200 in the lower right panel, for all subhaloes in 
the lower resolution ACDM simulation. The comoving trans- 
verse separation Ar± as a fraction of -R200 for the subhaloes 
that are selected by 9 < # max is shown as a red hashed re- 
gion while the distribution of those not selected is shown in 
blue. In the upper right panel in Fig. [3] the distribution of 
R200 for the parent haloes is plotted in green. Selecting pairs 
according to # max gives rise to a sample containing most of 
the subhalo pairs which have Ar± < R200 , with only a small 
number of pairs with Ar± > -R200 that happen to lie at an 
angle 9 < <9 max . 

We find the following selection rules give rise to a pop- 
ulation of pairs whose measured moment matches the pre- 
dictions of the AAP function extremely well, provided the 
correct cosmological model is assumed (see Section [5. ip : 

• The upper limit of the line of sight velocity difference 
should correspond to retaining 95% of the total distribution 
of pairs in the sample (grey shaded region in the left panel 
of FigHJ). 

• The maximum observed separation of a pair, 9 m „, 
should correspond to retaining 50 — 60% of the distribution 
for all pairs. 

Subhalo pairs in redshift space chosen according to the 
two selection criterion given above give rise to the blue 
hashed region shown in Fig. For this ACDM simulation at 
z = this corresponds to approximately 19,000 pairs with 
Av < 950 km/s (95% of the distribution) and 9 < 6.5 x 10" 4 
rad (50% of the distribution), bf Note these specific values 
quoted for Av and 9 are only for illustration. The selec- 
tion criteria presented in the two bullet points above should 
be applied the to parent sample of galaxy pairs when im- 
plementing this test. If we select pairs by restricting Av 
only, we retain 38,000 subhalo pairs and then the differ- 
ence between the measured mean and the corresponding 
AAP function is approximately 30%. This distribution is 
shown in green in Fig. [3] Selecting pairs with 9 < # max 
and no restriction on Av, i.e. including correlated pairs with 
(Au||/Ar) 7^ 0, results in a mean that differs from the cor- 
responding AAP function by less than 1%, provided the full 
expression for a in Eq. 1 is used (shown in orange in Fig. - 
Note if the full expression is not used then the difference is 
4%. We find that the measured moment is most sensitive to 
the first two selection criterion chosen by Marinoni & Buzzi 
and relatively insensitive to the minimum comoving sepa- 
ration of the pair and the comoving distance to the nearest 
neighbour. Note this is partly because we only consider pairs 
of subhalos from the same halo. 

In Fig. [5] we plot the measured mean of each of the 
distributions shown in Fig. [3] as a function of redshift. The 
red dot-dashed line shows the mean of the distribution of all 
subhalo pairs in redshift space with no cuts. By restricting 
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Figure 5. The sample means of the anisotropic distributions of 
pair tilt angles shown in Fig.[3]for subhaloes in the low resolution 
ACDM simulation, as a function of redshift. The mean of the dis- 
tribution of all pairs in redshift space (no selection cuts) is shown 
as a red dot-dashed line. The measured mean for pairs selected 
with only a cut in At) (8) is shown as a green dot-dot-dashed 
(orange dotted) line. Once we impose a cut in both Av and 9, 
the measured mean of the distribution (blue dashed line) agrees 
with the corresponding predicted AAP function (solid black line). 
The light green shaded region shows the uncertainty on this pre- 
diction because we have measured a from the simulation which 
has a finite number of pairs. The AAP function plotted here was 
found assuming a ACDM cosmology for H(z). The error bars on 
the data points are estimated from a jack knife sampling of the 
subhalo pairs using 100 subsamples of the data. 

the sample, using either a cut in Ad or 9, we obtain the mean 
shown as the green dot-dot-dashed line or orange dotted line 
respectively. Selecting subhalo pairs according to the two 
selection cuts discussed above results in a measured mean 
(blue dashed line) which is in very good agreement with the 
predictions of the corresponding AAP function (solid black 
line) when we measure a directly from the simulation. Note 
each distribution has its own associated AAP function, with 
a normalization set by the pairs in each sample. 

The error bars in Fig. [S] have been calculated by jack- 
knife sampling the subhalo pairs by grouping the data into 
100 sets containing equal numbers of subhaloes, and then 
successively removing one set at a time, calculating the sam- 
ple mean for the remaining haloes and computing the vari- 
ance amongst the measured means (see iNorberg et alll2009 , 
for a discussion of the reliability of the jackknife technique). 
We have verified that these errors change the AAP function 
by less 1% if we vary the sample size to 25 or 50 subsamples 
at a given redshift. The error on the AAP function, shown 
as a green shaded region in Fig. [5] was found using a similar 
method to find the variance in a at each redshift splitting 
the pairs in the simulations into 100 subsamples. The er- 
rors on both (sin 2 t) and on the AAP function increase with 
increasing redshift as the number of pairs decreases. This 
happens because of the fixed resolution of the simulation, 
which means that we resolve a progressively smaller frac- 



tion of the subhalo population with increasing redshift. A 
similar drop in the number of pairs would happen in a flux 
limited galaxy survey. 

We have tested the stability of the method by compar- 
ing simulations of different resolution. The sample mean, 
(sin 2 r) (Eq. [l|, of the distribution of subhalo pairs in red- 
shift space in the two ACDM simulations, higher and lower 
resolution, is shown in Fig. HJ] at different redshifts. The 
mean for the lower resolution simulation is shown as blue 
squares in Fig. [6] The AAP function using the measured 
value for a at each redshift is shown as a solid black line as 
in Fig. [5] The sample mean from the higher resolution simu- 
lation is shown as purple circles with the corresponding AAP 
function shown as a solid light blue line. In both the lower 
and higher resolution simulations, the measured distribution 
of pairs shows excellent agreement with the predictions of 
the AAP function assuming a ACDM cosmology, and also 
agree with each other in shape and amplitude within the 
error bars. The difference between the AAP function for the 
higher resolution (blue shading) and for the lower resolu- 
tion simulation (solid black line) in Fig. [6] is also due to the 
difference in resolution between the two simulations. If we 
select only subhalos from the higher resolution simulation 
that have a halo mass of M ^ 9 x 10 12 /i~ 1 Mq, which cor- 
responds to the minimum halo mass selected by the FOF 
algorithm in the lower resolution simulation, we obtain the 
red stars with errors bars plotted in Fig. [6] These points 
are almost coincident with the corresponding measurement 
from the lower resolution simulation (blue squares) , agreeing 
to better than la. 

We also make contact with an observational galaxy sam- 
ple in Fig. |6] If we select subhalos from main halos which 
have a mass of M ^ 1 x 10 14 /i _1 Mq we obtain the mean 
plotted as orange crosses in Fig. [6] Again these results are 
consistent with the means for the lower resolution simula- 
tion at each redshift. This mass corresponds to the min- 
imum halo mass expected to contai n two or more lum i- 
nous red galaxies (LRGs) on average (| Almeida et al.ll2008l ). 
This subhalo selection is relevant for a spectro scopic redshift 
surve y such as the SDSS-III BOSS survey (|Schlegel et al.l 
2007) which will target LRGs in the redshift range z < 0.7. 
Without applying any selection cuts, we find approximately 
27,000 subhalo pairs at z = which share a common 
halo of JO lx 10 14 /i -1 Mq, at z = 0.25 and z = 0.5 
the number of subhalo pairs are approximately 18,000 and 
11,000 respe c tively . From the first semester of BOSS data, 
IWhite et al.l (|201ll ) estimate that the cumulative probabil- 
ity that a galaxy in their sample is hosted by a halo of mass 
M > 1 x 10 14 /i _1 Mq is about 5%. If we extend this prob- 
ability to the full sample of LRGs expected by BOSS with 
space density ft — 3 x 10~ 4 /i 3 /Mpc 3 , this corresponds to 
approximately 13,000 pairs of LRGs in the redshift range 
z — 0.5 — 0.6. This is similar to the number of pairs we 
obtain from the higher resolution simulation restricting to 
halos with to M>lx 10 14 ft _1 MQ at z = 0.5, shown by 
the orange crosses in Fig. [6] 

The errors on the AAP function as measured by Mari- 
nom & Buzzi, a = 5.791J5; 3 , 2 , are substantially larger than 
ours due to the uncertainty in fitting for the parameter a at 
z = with a smaller number of pairs. Our higher resolution 
simulation has approximately 4 times more subhalo pairs 
than the lower resolution simulation, after making the selec- 
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Figure 6. The sample means of the anisotropic distribution 
of pairs of subhaloes measured in the low and high resolution 
ACDM simulations are shown as blue squares and purple circles 
respectively. The AAP function given in Eq. [9] using the value for 
a measured from the simulation at each redshift is shown as a light 
green and a blue shaded region for the low and high resolution 
simulation respectively. The mean measured using subhalos from 
the higher resolution simulation that have a parent halo mass of 
M>9x lQ la fc -1 Mg and M > 1 X W li h,- 1 M & are shown as red 
stars and orange crosses respectively. Note that for the measured 
mean of subhalos within a parent halo of M ^ 1 X 10 14 h~ 1 Mq 
at 2 = 2 there are only 49 pairs and we used 10 subsamples to 
find the jackknife errors. 



tion cuts discussed in Section 4, which gives rise to error bars 
which are approximately 50% smaller in the higher resolu- 
tion run (see Fig. [6]). The sample of pairs used by Marinoni 
& Buzzi is approximately 25 times smaller then the sample 
from our lower resolution simulation. 

We have verified that by applying the Marinoni & Buzzi 
selection cuts to our parent sample of subhalo pairs in the 
lower resolution simulation gives a = 5.69, which is consis- 
tent with the value for a obtained by these authors. How- 
ever, we find that the measured mean for this simulation 
sample does not agree with the AAP function within the 
error bars (if our sample was the same size as that used by 
Marinoni & Buzzi our errors would be significantly larger 
and the two would agree in this case). This demonstrates 
the need for the robust resolution independent selection cri- 
teria we have presented here. 



5 APPLICATION: A NEW TEST 

In this section we use the selection criteria outlined in Sec- 
tion 4 to test the predictions of the AAP function, Eq. [9] 
using the distribution of subhalo pair angles measured in 
N-body simulations. 

The accuracy of this test relies on two key variables: the 
cosmological expansion history assumed, H(z) and the nor- 
malization parameter, a — _r7 ~ 1 ((Ai; 2 )/(Ar 2 )) 1/ ' 2 . We con- 
sider the impact of uncertainties in each of these variables 
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Figure 7. The first moment of the anisotropic distribution of 
pairs of subhaloes measured in a ACDM and two quintessence 
dark energy simulations as a function of redshift. Measurements 
for ACDM, INV and SUGRA cosmologies are shown as blue 
squares, red circles and green triangles respectively. The AAP 
function using the measured value for a at each redshift for each 
cosmological model is shown as solid black, red and green lines for 
ACDM, SUGRA and INV respectively. The shaded bands show 
the uncertainty on the AAP function for each cosmology. Note 
that the error bars for the AAP function for SUGRA are similar 
to those for the INV model and are not shown for clarity. 



in turn. In Section [5T] we present the measured anisotropic 
distribution of the orientation of pairs, selected according to 
the prescription set out in Section 4, and its first moment at 
different redshifts together with the predicted distribution 
using the AAP function in a ACDM and in two quintessence 
dark energy cosmologies. In order to test the ability of the 
theoretical model to distinguish different cosmologies we will 
assume perfect knowledge of the correct H(z) and a in the 
first instance. We then consider how an observer would mea- 
sure a and the impact this has on the results, again assum- 
ing the correct H(z). We relax the assumptions further in 
Section I5.2I where an incorrect cosmological expansion his- 
tory is used to analyse the data. This is done by measuring 
the distribution of subhaloes in the INV and SUGRA dark 
energy simulations assuming a ACDM cosmology to infer 
distances to t he pair. We will show tha t the method, as im- 
plemented in iMarinoni fc Buzzil (|2010h . fails to exclude the 
wrong cosmology. Consequently, we propose a new method, 
which uses the theoretical model discussed so far but which 
exploits additional information about a from the numerical 
simulations. In Section 5.3 we show that this method can be 
successfully applied to test dark energy. 



5.1 Testing the method: theory versus 
observations 

First of all, we test how the approach discussed in Section 
2 can distinguish different cosmologies. We put ourselves in 
the idealised situation of an observer who knows the correct 
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cosmological model to compute distances and is able to mea- 
sure peculiar velocities precisely to find a at each redshift. 
In Fig. [7| the mean of the redshift space distributions of sub- 
halo pairs for ACDM and the two quintessence dark energy 
models are plotted as a function of redshift. The results for 
ACDM are the same as those shown in Fig. [15] for the lower 
resolution simulation. The measured sample mean for the 
INV dark energy model is shown as red-orange circles with 
error bars while the results for the SUGRA model are shown 
as green-grey triangles. The predicted AAP function for each 
of these models, using the correct expansion history and the 
value for a measured at each redshift, is shown as a solid red 
line for the INV model and a solid green line for the SUGRA 
model. The uncertainties on the AAP function are plotted 
as a red shaded region for the INV model. The errors for the 
SUGRA model are similar but are not plotted in Fig. [7] for 
clarity. The errors shown on both models for the measured 
mean and the AAP function were found using an identi- 
cal jackknife sampling method to that used for the ACDM 
result. It is clear from Fig. [7] that the measured mean for 
the three simulations agree with the respective AAP func- 
tion, provided the correct expansion history is known and 
that the parameter a can be determined at each z. As these 
results show, the measurements in a ACDM or a dynami- 
cal dark energy model agree very well with the predictions, 
if the correct cosmology is used to analyse the data. For 
the two quintessence models the deviations from ACDM are 
due to the different expansion histories (see Fig. [2} ■ This is 
a consistency check which confirms that the method works. 

In reality, in a galaxy survey, it is not possible to mea- 
sure the parameter a accurately at high redshifts because 
of the difficulties associated with measuring galaxy peculiar 
velocities to sufficient precision. We shall now degrade the 
status of the idealised observer mentioned above and con- 
sider a more realistic observer who still knows the correct 
cosmological model but who is unable to measure a directly 
at any redshift other than z = 0. Using the measured dis- 
tribution of pairs at z = 0, we fit the distribution given in 
Eq. [5] to set a and test the accuracy of the AAP function 
using this a(z = 0) value at each redshift, as suggested by 
Marinoni & Buzzi. If a does not evolve with redshift we 
would expect this approach to result in accurate agreement 
between the measured mean and the AAP function. 

In Fig.|8]the measured distribution of the angle r, in ra- 
dians, for ACDM is shown as a red hashed region with error 
bars. Note the y-axis shows the fraction of the total number 
of pairs per bin. The distribution (Eq. [6} with best fitting 
value for « F it = 5.67 ±0.1 (with lcr errors) is plotted as a 
purple dashed line. The grey dotted lines show the distribu- 
tion adopting a + la and a — la. Note the error we obtain 
for a, 0.1, is much smaller than that obtained by Marinoni 
& Buzzi (0.3) due to the difference in sample size and the 
different methods used to estimate the errors. This value for 
a agrees with the measured value from the simulations of 
a = 5.56. In Fig. [9] the AAP function assuming a ACDM 
cosmology and using this a F i T (z = 0) value at each redshift 
is shown as a black dashed line with error bars. The mis- 
match between this curve and the simulation results clearly 
indicates that a evolves with redshift, invalidating one of 
the main assumptions made in the analysis of Marinoni & 
Buzzi. Note that this black dashed line in Fig. [9] is much 
smoother than the shaded green band for the AAP function 
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Figure 8. The fractional distribution of the angle r in radians 
of pairs of subhaloes measured in ACDM at z = 0. The error bars 
on each bin are calculated by jackknife sampling after dividing 
the catalogue of subhalo pairs into 100 subsamples and calculat- 
ing the variance amongst the distributions measured after suc- 
cessively removing one subsample at a time. The purple dashed 
line shows the distribution in Eq. \E\ with the best fit value of the 
normalization parameter cifit = 5.67 ±0.1. The grey dotted lines 
show the 1 <r error on the best fit distribution. 

in ACDM where the value of a is measured directly from 
the simulation at each redshift. Using the z = value for 
a produces an AAP function which is systematically and 
significantly below the measured results for a ACDM cos- 
mology for z > 0. Applying this measure of a as proposed 
by Marinoni & Buzzi could lead to a spurious detection of 
deviations from ACDM. 

It is clear from Fig. [9] that the method proposed by 
Marinoni & Buzzi contains a serious systematic error which 
is apparent when applied to a large sample of pairs. Marinoni 
& Buzzi considered a smaller sample than in the simulations 
where the statistical errors dominated this systematic. 

It is clear from Fig. that a does evolve with redshift 
and that we can improve on the estimates of this parameter 
by fitting Eq.|6]to the measured distribution at each redshift. 
In Fig.[!J]the AAP function in a ACDM cosmology using the 
best fit values for a measured at each redshift is shown as a 
red dot-dashed line with error bars. The jackknife errors on 
a are estimated using 100 subsamples for the distributions 
at z — — 1 and 50 subsamples for z = 1.5 and 2 as there 
are fewer pairs at these higher redshifts. This approach to 
measuring a gives much better agreement with the mean 
measured from the simulations, shown as blue-yellow data 
points in Fig. Note this method of extracting a assumes 
that the correct cosmology is ACDM. 

5.2 The test assuming a particular cosmology 

In this section we are no longer idealised observers who know 
the correct cosmology, so the only possible choice is to as- 
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Figure 9. The first moment of the anisotropic distribution of 
pairs of subhaloes measured in a ACDM simulation (blue-yellow 
squares) as shown in Fig. [7] The AAP function using the 2 = 
best fit value for a at each redshift (see Fig. [SJ is shown as a 
dashed black line with error bars. If we fit for a using the mea- 
sured distribution at each redshift we obtain the red dashed line 
prediction for the AAP function assuming a ACDM cosmology. 
The error bars on the black dashed and red dot-dashed line are 
the 1 cr errors obtained by fitting for a. 



Figure 10. The measured mean of the anisotropic distribution 
of pairs of subhaloes in the INV dark energy simulation assum- 
ing a ACDM cosmology to find the comoving distance to each 
pair member (red circles). The purple dot-dashed line shows the 
predicted AAP function, assuming a ACDM cosmology, and us- 
ing the best fit a found at each redshift. The measured mean 
for the SUGRA model, assuming a ACDM cosmology, is shown 
as green triangles and the predicted AAP function, assuming a 
ACDM cosmology, and using the best fit a found at each redshift 
is shown as a grey dashed line. 



sume the same cosmology in the data fitting and in the theo- 
retical prediction of the galaxy distribution. Specifically we 
will assume ACDM, for simplicity, in order to set the ex- 
pansion history H(z) in Eq. 8 and to compute the comoving 
distances in Eq. [1] as well as to extract the parameter a. In 
order to find the parameter a we must fit to the observed dis- 
tribution of the orientations of pairs which has been found 
also by assuming a ACDM cosmology. Assuming that the 
true cosmological model chosen by nature is a dynamical 
dark energy model, for instance the INV or SUGRA cos- 
mology, we will check if the wrong cosmology, ACDM in our 
case, can be excluded or not, and consequently if the method 
is applicable to future galaxy surveys. For this analysis we 
take subhalo pairs in the INV and SUGRA simulations and 
at each redshift we rescale the comoving distances to match 
those which would be computed by an observer assuming a 
ACDM cosmology. 

The ensemble average of Eq. [T]for each subhalo pair in 
the INV (SUGRA) simulation is plotted in Fig. \W\ as red 
circles (green triangles), with error bars. If we fit for a at 
each redshift, we obtain the purple dot-dashed line in Fig. llQI 
for the INV model. Although we have assumed, incorrectly 
a ACDM cosmology we find that the AAP function agrees 
with the measured sample mean for the INV model at each 
redshift. In Fig. \W\ a similar analysis is presented for the 
SUGRA model. The measured mean for this dark energy 
model, assuming a ACDM cosmology to compute comoving 
distances, is shown as green triangles. The AAP function 
using the best fit value for a at each redshift and a ACDM 
expansion history is shown as a grey dashed line. Again, 
theory and observations agree when we would expect them 



not to as we have used the wrong cosmology in the AAP 
function and to compute distances. Our results show that 
the AAP function, using either a fixed value of a(z — 0) or a 
best fit value at each redshift, is not an accurate model with 
which to test for dynamical dark energy models if the correct 
cosmological model is unknown and that further input from 
numerical simulations is needed to arrive at a viable test. 

5.3 A blueprint for probing dark energy 

From the previous section it is clear that the cosmological 
test proposed by Marinoni & Buzzi relies heavily on mea- 
suring the parameter a accurately at each redshift. The as- 
sumption that a does not vary with redshift is incorrect 
and could falsely rule out ACDM if this test is misapplied 
to pairs of galaxies in future surveys. The value of a also 
depends on the cosmological model. For example at z — 
the values for the ACDM, SUGRA and INV cosmologies are 
a = 5.56, a — 6.32 and a — 5.32 respectively, with a typical 
error of 0.1. The difficulty is not just a problem of measur- 
ing a accurately but stems from the fact that the cosmology 
assumed affects both the data and the theoretical prediction 
in a way which cannot be disentangled. The accuracy and 
predictive power of the AAP function can be restored if in- 
stead of measuring a from observations, we employ N-body 
simulations which contain a comparable number of subhalo 
pairs to the number of galaxy pairs in the survey under con- 
sideration. It is clear that independent information about a 
is necessary and numerical simulations play an important 
role in providing these predictions in a given cosmology. We 
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propose a new approach to measuring dark energy, where 
observational measurements of the mean of the anisotropic 
distribution of pairs and predictions of the AAP function 
from numerical simulations are combined. 

The new method we propose to test a given cosmology 
is as follows: 

• An observer assumes the cosmology to obtain the co- 
moving distances needed to calculate the ensemble average 
of Eq. Q] for a sample of pairs of galaxies, selected using the 
criterion in Section 4, at different redshifts. 

• Using a N-body simulation of the same assumed cos- 
mology and with a comparable number density of pairs and 
volume to the galaxy survey, the observer can then construct 
a similar catalogue of pairs according to Section 4 and find 
the value of a at each redshift. 

• This gives rise to a prediction for the AAP function 
which can be compared with the means measured from the 
galaxy survey at each redshift and the assumed cosmology 
can be verified or excluded. 

Note if the AAP function measured from the simulation 
and the measured mean of the galaxy pair sample, analysed 
assuming the same cosmology disagree, then a suite of N- 
body simulations of different cosmologies would need to be 
run. The AAP function from each simulation should then 
be compared to the measured mean, computed assuming 
the cosmology used in the simulation. This test is realistic 
given current computing resources. 

In Fig. [11] we use the INV and SUGRA simulations to il- 
lustrate this method. In the upper and lower panels we show 
the measured means for the INV and the SUGRA dark en- 
ergy simulations respectively, which are treated here as the 
'observed' pair sample. In this example we are testing a 
ACDM cosmology and use it to compute the distances in 
each case, as in Fig. 1101 together with the predicted AAP 
function from an N-body simulation of ACDM where a is 
measured directly from the simulation (green shaded re- 
gion). It is clear that for z < 1 the INV model and the 
SUGRA model can be distinguished from the AAP function 
predicted in a ACDM cosmology. This result shows that if a 
SUGRA or INV model is the correct cosmology for the Uni- 
verse then ACDM can be ruled out. If there is a mismatch 
between the measurement from the observed pair sample 
and the simulation calibrated AAP prediction as in Fig. 1111 
then a new simulation with a different expansion history is 
required until an acceptable match is found. 



6 SUMMARY 

The distribution of the orientation of pairs of galaxies is uni- 
form in real space in a homogeneous and isotropic universe. 
However in redshift space, two effects lead to the inferred 
distribution becoming skewed. Firstly, an observer has to 
assume a cosmology to convert positions on the sky and 
redshifts into distances. A mismatch between the assumed 
and underlying cosmology introduces an error in the radial 
distance to a galaxy. Secondly, peculiar motions introduce 
distortions which break the connection between the mea- 
sured redshift and the actual distance. Both effects result in 
an apparent displacement of galaxies along the line of sight. 
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Figure 11. Upper panel: The measured mean of the anisotropic 
distribution of pairs of subhaloes in the INV dark energy simula- 
tion assuming a ACDM cosmology to find the comoving distance 
to each pair member (red circles) . The predicted AAP function for 
a ACDM cosmology, using the value of a measured directly from 
the lower resolution N-body simulation, is shown as a solid black 
line. Lower panel: same as the upper panel but for the SUGRA 
quintessence model assuming a ACDM cosmology to determine 
comoving distances (green triangles). 



iMarinoni fc Buzzil (|2010h proposed that the distortion 
of the distribution of the angle subtended between galaxy 
pairs as viewed in redshift space can be modelled by a sim- 
ple Doppler shift in the galaxy positions. This procedure 
gives rise to a theoretical prediction for the distribution in 
redshift space, which is referred to as the AAP function. 
We have tested the accuracy of this model using subhalo 
pairs identified in N-body simulations of cosmologies with 
different dark energy models. 

The AAP function depends on two variables: the "nor- 
malization" parameter a = Hq 1 ((Av 2 /Ar 2 )) 1 / 2 , and the 
expansion history, H(z), which depends on the cosmol- 
ogy. In this paper we present the AAP function normal- 
ized in three different ways: (i) using the relation a = 
Hq 1 ((Av^/ Ar 2 )) 1 / 2 , we can measure a directly from the 
simulation at each redshift, (ii) we can measure a at z = by 
fitting to the measured distribution and then, assume that 
a does not evolve with redshift, and apply the z = nor- 
malization to specify the mean of the distribution of pairs at 
different redshifts (as suggested by Marinoni & Buzzi), (iii) 
we can apply case (ii) but fitting for a at each redshift using 
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the measured distribution and not just at z = 0. When we 
measure a directly (case i), we obtain excellent agreement 
between measurements of the mean from the simulation and 
the predictions of the AAP function. If instead we retain the 
best fit z = value, a(z = 0), at each redshift (case ii), we 
do not find a good match between the theory and the simu- 
lation measurements. This demonstrates that simply fitting 
for a at z — and assuming it does not evolve with redshift 
is not accurate. In fact such an approach would incorrectly 
rule out the cosmology used in the simulation. If we fit for 
a at each redshift using the simulation (case iii), then we 
again recover an excellent match between the theory and 
simulation results. 

We use a large sample of subhalo pairs which do not 
necessarily reside in fully relaxed and virialised haloes that 
have detached from the Hubble flow. This is demonstrated 
by the fact that we measure a different value for a in different 
cosmologies, see Section 5. 

Note that each of the above cases consider idealised ob- 
servers who know the correct cosmology to compute dis- 
tances and H(z). The measured mean of the distribution of 
pair angles (Eq. [TJ depends on the cosmological model as- 
sumed to convert position on the sky and redshift to comov- 
ing distance. The AAP function also depends on cosmology 
through H(z). As a result, the measured mean and the AAP 
function will not agree if the wrong cosmology is assumed 
(the Alcock-Paczynski effect). Using two quintessence dark 
energy simulations (labelled INV and SUGRA), we have 
tested if the AAP function reproduces the measured mean 
of the distribution when, in the first instance, we know the 
correct cosmology (the 'perfect' observer case), and in the 
second instance, when we instead assume ACDM (i.e. the 
'real' observer who has no prior knowledge of the under- 
lying cosmology). The two dark energy models we consider 
have an evolving equation of state which is compatible with 
current observations of the CMB, BAO and Type la SN dis- 
tances. We find that, for a perfect observer who knows H(z) 
and a exactly, the AAP function and the measured means 
are in very good agreement for both the SUGRA and the 
INV models. 

Consider now performing the same exercise using the 
SUGRA and INV simulations, as a real observer who does 
not know the underlying cosmology and so assumes a ACDM 
cosmology, and who uses the best fit value for a at each 
redshift. We might expect that the theory should not match 
the measured mean for the dynamical dark energy models. 
However, we find that, by fitting for a using the observed 
distribution in the simulations, we instead recover a model 
which incorrectly matches the observations extremely well 
for both dark energy cosmologies, even though we have as- 
sumed a ACDM model. The consequences are that, in a uni- 
verse with evolving dark energy, we would find that a ACDM 
model incorrectly matches the observations invalidating the 
methodology. 

In this paper we have proposed a new formulation of 
the test of Marinoni & Buzzi in which the distribution of 
galaxy pairs can be analysed without prior knowledge of the 
cosmology. The measured distribution of angles should be 
compared with predictions for the AAP function using a 
reference N-body simulation to directly measure a. We have 
shown that the subhalo pairs in two quintessence dark en- 
ergy simulations, which are treated as the "observed" pair 



sample in this instance, produce a different measured distri- 
bution to that predicted in a ACDM simulation even when 
analysed after assuming (incorrectly) a ACDM cosmology. 
In the new test, the AAP function is normalized with ref- 
erence to a simulation with the same cosmology assumed 
to analyse the observations. The predicted AAP function 
and measurement will only agree if the assumed cosmology 
matches the true cosmology. If this is not the case, then a 
new reference simulation must be generated with a revised 
expansion history, to see if an improved match to the ob- 
served distribution of galaxy pair angles can be obtained. 
We find that, by measuring the mean of the distribution as 
a function of redshift, we should be able to detect devia- 
tions from a ACDM expansion history at the level of 2% in 
a box of volume ~ 3/i~ 3 Gpc 3 . This new test complements 
the constraints on the present value of Hubble's parameter 
provided by observations of Type la SN which constrain Ho 
to ~ 3% l|Riess et al.ll201lh and improves on constraints of 
H (z) at higher redshifts which are accurate to about ~ 10%. 



7 CONCLUSIONS 

Distinguishing between competing scenarios for the acceler- 
ating expansion of the universe is a major challenge for both 
observational and theoretical cosmologists. The expansion 
history and distance-redshift relations are remarkably close 
between viable models which satisfy the currently available 
constraints. A convincing determination of the nature of 
dark en ergy will require a c ombination of probes for two 
reasons (jAlbrecht et al.ll2006h . Firstly, the small differences 
in the expected signals from a given probe mean that sys- 
tematic effects become important. Applying different probes 
will allow us to see whether or not a measured signal is ro- 
bust to systematics. Secondly, some existing tests cannot 
distinguish between some classes of dark energy model. New 
probes are therefore needed to break such degeneracies. 

We have tested one such example of a new probe, the 
distribution of angles subtended between pairs of galaxies. 
This distribution is distorted by the peculiar motions of 
galaxies and also by the choice of cosmology adopted to 
transform observed positions into comoving distances. The 
origins of this test can arguably be trace d back to Alcoc k 
fc Paczynski (1979 ), and it was refined bv lPhillippsI il 19941) . 
iMarinoni fc Buzzll (|2010l ) applied the test to the angle be- 
tween pairs of galaxies and crucially included redshift space 
distortions. 

We have used numerical simulations of structure for- 
mation to assess the performance of the test. The mean 
of the distribution of pair angles varies with redshift and, 
furthermore, is measurably different between cosmologies. 
A comparison between a theoretical model for the pair an- 
gle distribution and the measurements from the simulations 
shows that the test, as originally proposed, is limited. The 
theoretical calculation requires a parameter to be specified 
to normalize the distribution of pair angles. Our simulations 
show that this parameter is redshift and cosmology depen- 
dent. 

It is possible to estimate the normalization of the pair 
angle distribution observationally, at redshifts z > 0, if the 
peculiar velocities of galaxies can be measured. For exam- 
ple, it was recently argued that accurate mean pairwise ve- 
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locities of pairs of Type la SN can be obtained by com- 
bining photometry from a survey such as Pan-STARRS 
Kaiser et al.ll2010 | ) or the Large Syn optic Survey Telescope 



LSST Science C ollaborations 2009) with follow-up spec- 



troscopy 1 Bhattacharva et alj|2011f ). At present the accu- 
racy of measurements of the peculiar velocity field is not 
adequate to distinguish between the models compared in 
this paper. 

Our proposed methodology avoids this problem by us- 
ing an N-body simulation with a similar number of pairs 
to the observational sample to normalize the distribution of 
angles. This secures the crucial step of setting the normal- 
ization of the theoretical distribution at each redshift. The 
detailed selection of the N-body sample of subhalo pairs is 
not important, avoiding the need to combine the simula- 
tion with a galaxy formation model. Furthermore, we have 
demonstrated that it is not necessary to have a knowledge of 
the true underlying background cosmology for the successful 
application of the test. 

The new method we have proposed is a powerful com- 
plement and extension to existing probes of dark energy. 
This is demonstrated by the ability of the pair distribu- 
tion to distinguish between cosmologies that cannot be sep- 
arated through the appearance of BAO or through the halo 
mass function. The technique can be applied already to on- 
going surveys, such the SDSS-BOSS survey l|Schlegel et al.l 
120071 ). and should yield competitive constraints. The method 
should also produce distinct signals for dark energy and 
modified gravity models which have identical expansion his- 
tories, through the different peculiar motions induced. 
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